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(N ■ ABSTRACT 

pL| I We have discovered an axially symmetric, well-defined shell of material in the 

^ I constellation of Cepheus, based on imaging acquired as part of the Galactic First 

Look Survey with the Spitzer Space Telescope. The 86" x 75" object exhibits 



> 



brightened limbs on the minor axis, and is clearly visible at 24 fim, but is not 



CN I detected in the 3.6, 4.5, 5.8, 8.0, 70, or 160 /xm images. FoUowup with 7.5 - 

40 /xm spectroscopy reveals the shell to be composed entirely of ionized gas, and 
CN I that the 24 fim imaging traces solely [O iv] 25.89/im emission. The spectrum also 

\^ [ exhibits weaker [Ne ill], [S ill], and very weak [Ne v] emission. No emission from 

^ I warm dust is detected. Spectral cuts through the center of the shell and at the 

'-^ ' northern limb are highly consistent with each other. The progenitor is not readily 

^ I identified, but with scaling arguments and comparison to well-known examples 

^ I of evolved stellar objects, we find the observations to be most straightforward to 

^ I interpret in terms of a young supernova remnant located at a distance of at least 

10 kpc, some 400 pc above the Galactic disk. If confirmed, this would be the 
first SNR discovered initially at infrared wavelengths. 



Subject headings: Subject headings: ISM: individual (SSTGFLS J222557+601148) 
— infrared: ISM — (stars:) circumstellar matter 



1. Introduction 

Immediately following the science verification phase of the Spitzer Space Telescope 
(Werner et al. 2004), a Galactic First Look Survey (GFLS) was performed as an early 
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demonstration of Spitzer's capability to do Galactic science with mid-far-IR imaging de- 
vices. The Infrared Array Camera (IRAC, described by Fazio et al. 2004), with channels 
at 3.6, 4.5, 5.8, and 8.0 /xm, and the Multiband Imaging Photometer for Spitzer (MIPS, 
described by Rieke et al. 2004), with channels at 24, 70, and 160 /xm. Observations from the 
GFLS^ revealed to us a shell-like object in Cepheus, roughly circular in geometry, 86" x 75", 
and possibly stellar in origin with an infrared morphology suggestive of a planetary nebula 
(PNe) or supernova remnant (SNR). However, the object (SSTGFLS J222557+601148) has 
unusual infrared colors. It is easily detected in the MIPS 24/im imaging, but is undetected 
in the other 6 imaging bands. Exhaustive searches through the literature and ground- and 
space-based facility archives did not provide any information about this source. No pointed 
observations have been recorded, nor is the shell detected in the 2 Micron All Sky Survey 
(2MASS), the (optical) Digitized Sky Survey (DSS), the Canadian Galactic Plane Survey 
(at 1420 and 480 MHz), or the ROSAT All Sky Survey. This region is not covered in the 
archives of the Chandra X-ray Observatory, the XMM-Newton Observatory, the Einstein 
Observatory, or the Very Large Array. Soon after completion of the GELS, spectroscopy 
covering 7.5 - 40 //m at low spectral resolution {R = A/AA ~ 75 — 125) was obtained with 
Spitzer's Infrared Spectrograph (IRS; described by Houck et al. 2004). The resulting IRS 
spectroscopy provides an immediate explanation for the unusual IRAC- MIPS colors, but not 
unambiguously for the nature of the shell or its progenitor. 

In this Letter we present the Spitzer ima^ng and spectroscopy of the shell source, point- 
ing out its unique Spitzer-selected mid-infrared properties. The observations described here 
provide a view of an elusive, if not unusual, category of evolved stellar object and a valuable 
observational perspective on the interpretation of broad-band imaging in conjunction with 
spectroscopy. 

2. Spitzer Observations and Results 

The IRAC and MIPS observations were obtained on 2003 Dec 8 (AOR keys 4959488 
and 4961280), and processed in the Sll.0.2 pipelines at the Spitzer Science Center. The total 
exposure time on the shell averaged 48 seconds for IRAC (comprised of 4 dithers of 12 sec high 
dynamic range exposures). The majority of the GELS had an exposure time of 60 seconds 
for the 24 and 70 fim channels, but the shell was close to the edge of the area surveyed in fast 
scan mode and instead had an average exposure time of 36 sec at 24 fim and 15 sec at 70 fira. 
Only the western ~50% (roughly) of the shell visible at 24 /xm was imaged at 70 and 160 /xm. 



^General information about the Spitzer FhS can be viewed at http://ssc.spitzer.caltech.edu/fls. 
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The l-(7 noise levels in the mosaics are < 0.08 MJy sr~^ in the IRAC bands, and 0.1, 2.0, 

and 6.0 MJy sr~^ at 24, 70, and 160 fim, respectively. The IRS low resolution spectra were 
obtained in Staring Mode on 2004 July 15 (AOR key 10066176) at two positions: through 
the center of the shell, and on the brightened northern limb, using the SL (only in the first 
spectral order covering 7.5 - 14.5 /xm) and LL (14.5 - 37.5 /im) modules with total exposure 
times of 240 s and 360 s, respectively. Background corrections for the LL spectra were 
performed using alternated off-source subshts, while the SL spectra were corrected with off- 
source measurements taken ~ 4'.4 to the north (AOR key 10065408) using the same on-source 
SL integration time. The basic cahbrated 2-D spectrograms were produced using the Sll.0.2 
pipeline. The spectra were extracted and photometrically calibrated offline using extended 
source calibrations derived from the zodiacal emission and a scheme of pixel weighting over 
the illumination profile of the extraction apertures. Absolute photometric uncertainties are 
estimated to be 15 - 20%, tied closely to the uncertainties in the calibrated COBE/DIRBE- 
based models of the zodiacal emission^ at IRS wavelengths. The relative spectrophotometric 
quality of these data is quite good, estimated to provide equivalent point source detections 
of 15 mJy at 17//m and 40 mJy at 33 /xm (3 a), so that all labelled features are salient. 
However, background corrections have left slight residuals of the diffuse H2 S(l) 17.0 //m 
and S(0) 28.3 iim because of an E-W gradient in background levels through the wavelength 
subslits of LL2 and LLl. 

The 8.0 and 24 //m imaging is shown in Figure la,b. The 24-/xm shell exhibits bilateral 
symmetry and brightened limbs. This symmetry could be explained by several phenomena, 
including: (i) a bipolar outflow from a rotating progenitor, with the projected axis of rotation 
along the major axis; (ii) enhanced density or gravitational confinement of material in the 
plane of an orbiting companion projected along the minor axis; and [iii) confinement of 
ejected material from the acceleration of energetic particles along the minor axis oriented 
perpendicular to the interstellar magnetic field. The northern limb appears brighter than 
the southern limb, but this could be due to inhomogeneities in the ISM dust emission along 
the fine of sight. Arcs of emission at 24 //m may extend to the east and west around the 
shell, possibly as an outer shell of previously expelled or swept up material. However, the 
relative orientation of these features is not consistent with that of the brightened limbs, and 
the arcs may be unrelated structures in the ISM. 

The distance to the shell is unknown (see Sec. 3.2). There is no evidence from the ISM 
morphology that the shell is interacting with the nearby carbon star V384 Cep (situated 
8'. 8 to the north of the shell) or any other nearby object in projection. The shell lacks an 



^http://ssc. spitzer.caltech.edu/docunients/background/node3. html. 
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obvious central source to identify with the progenitor or its remnant. Because of the lack of 
emission at 70 /xm or a strong 24 fim point-like source, we must exclude that the source of 
the shell is obscured by a disk or thick envelope of cool dust. 

In Figure 2 we present the 7.5-40 /im spectra taken through the center and northern 
limb of the shell, after subtraction of the diffuse background emission. Both positions exhibit 
simple and near-identical spectra, dominated by [O iv] (25.89/im), accompanied by low 
excitation lines of [Ne ill] (15.56 /im), [Ne v] (14.32, 24.32/im), and [S ill] (18.71, 33.48 fim). 
Line identifications and fluxes are listed in Table 1. The IRS spectra provide an immediate 
explanation for the IRAC+MIPS colors: no continuum or dust features are detected, so the 
24 yum emission arises almost exclusively from the [O iv] emission line, while non-detection 
in the IRAC 8-//m channel is corroborated (at IRS sensitivities) by the lack of any emission 
over the passband. Non-detections at 3.6, 4.5, and 5.8 /im sets a 1-a upper limit flux of < 
0.1 MJy sr~^ in these bands. 

The IRS spectra show little variation in ionic species or relative intensities of the hues 
between the center and northern limb. This indicates little difference in the density of 
material sampled at these positions with the IRS slits, and that the object is indeed only a 
shell of gas that cannot be discerned to extend radially inward. The measured line emission 
from the limb compared to the center is systematically weaker due to the smaller volume 
of material integrated in the spectral extraction. Using the ratio of the [S III] (33.5/18.7) 
emission line intensities, we estimate an electron density in the range of 2.5-4.0 (±1.0) x 10^ 
cm~^ for collisionally excited lines (Rubin et al. 2001). The higher value corresponds to 
the [S III] line ratio at the northern limb, but both values are within the measurement 
uncertainties. 

Aside from the rather restricted range of ionic species and low energy transitions, we 
are surprised to flnd no indications of thermal emission from dust grains in the gas shell at 
any IRS wavelength or at 70 /xm (with dust temperatures in the corresponding range of 50 
- 500 K). Dust has either been destroyed (through sputtering, for example), or could not be 
condensed in the outflow. We suggest that very little interaction between the material in 
the shell and the ISM has occurred, and that the shell has a very high gas/dust ratio, with 
a surface brightness dominated by [O iv] emission. 
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Fig. 1.— Spitzer GFLS mosaics covering SSTGFLS J222557+601148 at (a) 8.0 fim and (6) 
24 /xm, at angular resolutions of 2" and 6" respectively. The shell-type SMC SNR IE 0102.2- 
7219 (centered at a = 01h04m02s, 6 = -72°01'50" J2000.0) is also detected with Spitzer 
cameras only at 24 /im (see Sec. 3.2), shown (c) in this mosaic from the Spitzer archive 
(Stanimirovic et al. 2005). 



Table 1: Principal IRS emission lines'^. 
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Units of A and line intensities are /im and 10 ^'^ W cm ^. 
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Fig. 2. — IRS spectra of the shell through the center (upper, offset by +0.15 Jy) and 
northern limb (lower). The spectra have not been rectified with basehne fits, nor corrected 
for extinction. Dashed vertical lines indicate where the SLl, LL2, and LLl spectral fragments 
are joined. Weak residuals of H2 S(l) and S(0) emission following correction for the diffuse 
background are indicated by 'x' near 17 and 28 //m. 
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3. Discussion 

3.1. An AGB or Post-AGB Shell, or Evolved Massive Star Nebula? 

Spherical shells of circumstellar material are typical among PNe, but the colors and 
associated emission spectrum of the GFLS shell are quite atypical. First, only pure gas 
is detected in the shell, and indicated to be limited in radial extent, or density bound, as 
evidenced by the similar values of Nf. derived from the spectral cuts through center and 
northern limb. In contrast, the gas and dust in AGB shells and PNe are radially distributed 
with large variations in Nf. within the individual nebulae. This is expected of slow, protracted 
outflows, often occurring in multiple events with variations in composition (from C- to O- 
rich) of material in the resulting multi-layered shell (e.g.. Waters et al. 1998). The shell 
also lacks an infrared-bright central source and spectral features that characterize AGB 
shells and PNe, which are heated O- and C-bearing dust and often simple metal oxides 
condensed in a dense outflow, creating as strong thermal continuum over a range of grain 
temperatures. The infrared spectra of PNe (e.g., Beintema et al. 1996) also typically exhibit 
simple and complex H-bearing molecules or hydrocarbons and higher ionized species with 
more upper-lying transitions of [S iv], [Mg iii-v], and [Ar iii-v], which are diagnostic of 
X-ray emission. Transitions from [Si II] and [Ne II], in ground state and excited levels, are 
observed in prototype PNe such as NGCs 7027, 6543, and 6302 (e.g., Beintema et al. 1996; 
Bernard-Salas et al. 2003; Bernard-Salas & Tielens 2005), and in eUiptical PNe such as 
IC 418, IC 2165, and NGC 5882 (Pottasch et al. 2004). Their absence in the spectra of the 
GFLS shell excludes coUisional excitation in slow-moving shocks (< 50 km s~^). 

In most hot ionizing source can be identifled within recently-formed PNe, but 

this is not a good criterion for very old PNe with progenitors which have evolved to the white 
dwarf phase. If the progenitor is now a white dwarf, it would probably be fainter than 18th 
magnitude in V (see the McCook & Sion 1999 catalog of compiled white dwarf photometric 
observations), escaping detection in the DSS. For surface effective temperatures spanning 
the range ~3400 K - 2 xlO^ K, we would not expect to directly detect a white dwarf fainter 
than about V = 19 mag in the Spitzer imaging bands. Spectroscopically, the nebular gas 
should contain the hydrogen lifted away from the surface of the star in its AGB or post-AGB 
phases, prior to the star emerging onto the H-deflcient white dwarf cooling track, and this 
is not consistent with the H-free IRS spectrum. 

Similarly, in the absence of dust or characteristic high ionization nebular and stellar 
wind emission lines indicative of mild to advanced CNO-processed material, the object is 
not likely to be the nebula around an evolved massive star counterpart, such as an OB-type 
supergiant or related luminous blue variable (LBV) or Wolf-Rayet star. Dust is generally 
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present in circumstellar nebulae of LBVs (e.g., Voors et al. 2000; Morris et al. 1999), but if 
the central star is optically shrouded then the 24 /xm emission from the central region should 
be very strong. Otherwise, the infrared stellar emission line spectrum is easily recognized 
in any of these phases by broad lines dominated by H, He, C, O, and various fine structure 
lines formed at velocities up to several 10^ km s~^ (e.g., Morris et al. 2004). 

3.2. Remnant of a Type I supernova? 

The lack of previous detection in X-ray and radio surveys does not exclude the shell 
from formation by material expelled in a supernova explosion, if the object is relatively 
distant, and/or little interaction between the shell and the ISM has occurred. The angular 
size of the shell is reasonable for a Galactic SNR if we constrain its distance to an upper 
value of 10 kpc and its time-averaged expansion velocity v^xp to less than 3000 km/sec (since 
the emission lines are unresolved). Adopting v^xp — 2000 km/s, then the radius of the 
remnant r = 2.0 x (age/lO^yr) pc, and the distance D — 10.0 x (age/lO^r) kpc. Reducing 
the expansion velocity decreases the distance and age, which would be more difficult to 
reconcile with lack of detection. An outside distance of ~10 kpc, on the other hand, is more 
reasonable when considering the shell's relatively unperturbed appearance and the inferred 
lack of interaction with the ISM that could be explained by location at a large distance 
above the Galactic plane. Since the object is at a latitude of 6 = 2°. 25, it would be 400 pc 
above the Galactic plane at a distance of 10 kpc. If the object is a SNR, then it is among 
the youngest (in the same kinematic age group as IE 0102.2-7219 and Cas A), and one of 
the smallest. Only the SNR Gl. 9-1-0. 3 (also a shell-type remnant) is comparably sized, at a 
diameter of 1'.2, and it is also one of the faintest (0.6 Jy at 1 GHz) among the 231 Galactic 
SNRs cataloged by Green (2004). 

A compelling but cautious comparison to the bilaterally symmetric shell-type remnant 
SN1006 also shows remarkably similar morphologies (see especially the ROSAT HRI image 
in Figure 1 by Dubner et al. 2002), with emission properties at radio and X-ray wavelengths 
that would make detection difficult were it placed beyond its distance of 2 kpc (Dubner 
et al. 2002). The thermal shell of SN1006 exhibits two non-thermal limb-brightened arcs 
to the NW and SE (Willingale et al. 1996; Winkler & Long 1997; Dyer et al. 2001), 
expanding into a smooth low density local ISM according to H I observations by Dubner et 
al. (2002). The brightened limbs have been proposed to be the result of the acceleration 
of relativistic particles along the NW-SE axis of symmetry, oriented perpendicular to the 
interstellar magnetic field (Roger et al. 1988). The flux density at 843 MHz is 17.5 Jy (Roger 
et al. 1988), or ~0.7 Jy were it located at a distance of 10 kpc. No mid- infrared observations 
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of SN1006 have been performed to our knowledge. 

Concerning the infrared colors of the GFLS shell, Stanimirovic et al. (2005) have recently 
pubhshed 5pzteer imaging (using all seven imaging bands) of the young (~1000 yr), 0-rich 
SNR IE 0102.2-7219 (hereafter E0102) in the Small Magellanic Cloud (SMC), and hkewise 
report IR detection of the spherical shell remnant only in the 24 //m channel (see Fig. 1). 
Although Stanimirovic et al. did not acquire IRS spectroscopy, they attribute the emission 
at 24 ^m to a combination of [O iv] line emission and dust heated to 120 - 130 K, where the 
latter comprises at least 40% of the total emission and coincides with regions of forward and 
reverse-shocked material traced by Chandra X-ray observations in the 0.3 - 10 keV range. 
Our IRS spectroscopy of the GFLS shell reveals only [O iv] line emission in the MIPS 24 /xm 
passband. Dust may be destroyed in a SNR cavity by sputtering on timescales Tgput ~ a/n 
(yr), where a is the grain size in /im and n is the ambient gas density in cm~^ (Dwek & 
Werner 1981). Grain sizes in the range of 0.01 - 0.1 //m will be eroded in less than 100 years 
at relatively high gas densities n ~ 10^ cm~^, if the observed [S ill] (33/19) line ratio is 
accurately representative. 

The IRS line spectrum of the GFLS shell is very similar to that of the 0-rich SNR 
Cas A (Arendt et al. 1999), which exhibits the same elements and ionization levels, but 
also [Ne ll], [S iv], [Ar ii-iii], and [Si ll] variously in different regions and high velocity 
knots associated with ejected material from the supernova explosion. Shock velocities in 
the range of 150 - 200 km/s are consistent with the Cas A observations, whereas higher 
shock velocities are needed in the GFLS shell to inhibit populating low ionization states to 
their critical densities. Using plane parallel shock models generated by the MAPPINGSII 
code (Sutherland Sz Dopita 1993; Dopita & Sutherland 1995) (at solar abundances), we 
find that [Ne ll] becomes weak at t'shock > 450 km s~^, but is not completely eliminated 
(to below the IRS detection limits) before [Ne v] becomes strong as velocities exceed 510 
km s"-*^. Nonetheless, shock speeds in the range of 450 - 500 km s~^ produce qualitatively 
consistent line ratios. Fast shocks also lead to the destruction of small (0.01 - 0.02 /im) dust 
grains and the mantles of larger grains (McKee 1989). If the GFLS shell is a SNR, then 
the progenitor could likewise be a massive star characterized in its last stage of evolution by 
high O abundance and H deficiency, terminating as a Type Ib/c supernova. 

We cannot exclude that it is the remnant of a Type la explosion as a result of core 
instabilities set up by high accretion of mass onto a white dwarf from a companion that 
filled its Roche lobe, or the defiagration of a white dwarf or degenerate core of a moderate 
mass star. The hydrogen- deficient state of white dwarfs and the morphology of Type la 
remnants would be consistent with our observations in this scenario. Indeed, the remnant 
SN1006 is widely believed to have been a Type la event (e.g., Wu et al. 1983). However, 
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models predict that the bulk of the ejected material in Type la explosions is composed of 
iron, following decay of ^^Ni by electron capture to form ^^Fe, given enough time to interact 
with the ambient medium (Wu et al. 1983, 1997). 

Observations at radio wavelengths, and further Spitzer imaging and spectra are needed 
to further test the SNR scenario for non-thermal emission and additional dynamic and 
kinematic constraints. If confirmed, this object would be the first SNR discovered by its 
infrared properties. 

We appreciate the helpful comments from Dr.s J. Bernard-Salas, W. Reach, and our 
anonymous referee. This research is based (in part) on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of 
Technology under NASA contract 1407. 
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